Heat Transfer Operators Associated with Quantum Operations 

Q. Aksak and S. Turgut 
Department of Physics, 
Middle East Technical University, 
TR-06800, Ankara, Turkey 

Abstract 

Any quantum operation applied on a physical system is performed as a unitary transformation 
on a larger extended system. If the extension used is a heat bath in thermal equilibrium, the 
concomitant change in the state of the bath necessarily implies a heat exchange with it. The 
dependence of the average heat transferred to the bath on the initial state of the system can then 
be found from the expectation value of a hermitian operator, which is named as the heat transfer 
operator (HTO). The purpose of this article is the investigation of the relation between the HTOs 
and the associated quantum operations. Since, any given quantum operation on a system can 
be realized by different baths and unitaries, many different HTOs are possible for each quantum 
operation. On the other hand, there are also strong restrictions on the HTOs which arise from the 
unitarity of the transformations. The most important of these is the Landauer erasure principle. 
This article is concerned with the question of finding a complete set of restrictions on the HTOs that 
are associated with a given quantum operation. An answer to this question has been found only for 
a subset of quantum operations. For erasure operations, these characterizations are equivalent to 
the generalized Landauer erasure principle. For the case of generic quantum operations however, 
it appears that the HTOs obey further restrictions which cannot be obtained from the entropic 
restrictions of the generalized Landauer erasure principle. 
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I. INTRODUCTION 



A fundamental result in the thermodynamics of computation is Landauer's erasure prin- 
ciple (LEP), which places a lower bound on the average heat dumped to a bath when the 
information in a memory device is erased. The principle holds for processes that satisfy 
two essential features: The process is carried out in the same way independent of the in- 
formation stored in the device and it restores the device to a known standard state at the 
end. These processes are usually known as the Landauer erasure processes. Landauer has 
shown that if the device holds / bits of information, then during any Landauer erasure, the 
average amount of heat that must be dumped to the bath is at least IksT \n2 where ks is 
Boltzmann's constant and T is the temperature of the bath 

The principle lies at the heart of the resolution of Maxwell's demon paradox j2-4|, an 
important factor that stimulates continued interest on the thermodynamics of information 
processing. The demon has also been investigated quantum mechanically j4-7|. Also, several 



alternative proofs |8l-ll 1| and generalizations 
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14j | of LEP have been given and different 
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lej. The 



approaches aimed at quantum information processing have been developed 
crucial ingredient in all alternative approaches to LEP is the description of the operation 
on the device as a unitary transformation on a larger combined system of the device and 
the bath where initially the two are uncorrelated and the bath is in thermal equilibrium {23]. 
Unitarity of this transformation is an absolute necessity for all physical operations |l3] while 
the noncorrelation and equilibrium assumptions are well-suited for all possible experiments. 
Although this viewpoint is contested [isj, LEP directly follows from the unitarity: in a 
Landauer erasure, the transformation maps different initial device-states to the same final 
device-state. Since the unitary transformation has to be one-to-one, this can only happen 
if the state of the bath is also modified. Finally, as the bath was initially in equilibrium, 
changes in its state imply an increase in its energy. Exact calculations then show that the 
average increase in the energy of the bath must exceed the Landauer bound. 

Can there be further implications of the unitarity of the transformation on the energy 
exchange with the bath? A recent article gives an affirmative answer to this question for 
the case of processes that manipulate classical information 1^. That article presents an 
alternative formulation of LEP in terms of the detailed heat transfers into the bath. In this 
non-entropic approach, one takes the individual heats dumped to the bath for each initial 
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logical state as the main quantities of interest and places restrictions between them. For 
example, if qi is the heat dumped to the bath during a Landauer erasure when the initial 
logical state is i, then it can be shown that 



an inequality which is first derived by Szilard [lO]. The main advantage of this formulation 
is its independence of the coding of the information into the logical states of the device. The 
standard statement of LEP, however, is dependent on the coding: if the device is used for 
storing information in such a way that the logical state i is used with frequency pi, then the 
device has information storage capacity of / = — ^ • pi log2 Pi bits and the Landauer bound 
states that 



where q = YliPiQi is the average heat dumped to the bath when the information in the 
device is erased. 

Note that Eq. ([2]) is valid for any distribution p, hence it provides an infinite number 
of restrictions on q^. It can be shown that Eq. ([2]) holds for all possible distributions p 
if and only if Eq. ([T]) is satisfied. Thus, Eqs. (IT]) and (E]) are two equivalent formulations 
of LEP. The standard form of LEP, Eq. (JJj) is more profound as it directly connects two 
disparate disciplines, thermodynamics and information theory. Yet, Eq. ([1]) is also somewhat 
attractive as it is a concise inequality that captures all of the infinitely many inequalities in 
Eq. ([2]). Both inequalities can be generalized to arbitrary nondeterministic processes. It is 
observed that the generalization of Eq. ([1]) contains further restrictions on heat exchanges 
that cannot be derived from the entropic restrictions of the generalized LEP jl^ . 

The purpose of this article is to adapt that coding-independent approach to quantum 
operations applied on devices holding quantum information. The major difference between 
the classical and quantum cases is in the description of the state: instead of a discrete index 
i, we now have a density matrix that describes the state of the device. Hence, in place of the 
numbers qi, we now have a hermitian operator Q whose expectation value with the density 
matrix gives the average heat dumped to the bath. This operator will be called as the heat 
transfer operator (HTO). 

The unitary transformation can be thought as moving information around in the com- 
bined system of the device and the bath. The resulting information exchange between the 
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q > IksT In 2 
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device and the bath changes the state and hence the energy of the bath. The HTO essen- 
tially provides a measure of this information exchange in terms of the energy given to the 
bath. The nature of the operation applied on the device (e.g., whether the information in 
the device is erased completely or only partially, or whether randomness from the bath is 
introduced to the device or not) has a strong effect on the HTO. This article investigates 
that relation between quantum operations and the associated HTOs. This relation essen- 
tially follows from the unitarity of the transformation on the combined system. However, it 
appears that the unitarity property is unnecessarily too restrictive; all of the results in this 
article can also be obtained from a weaker assumption, namely that the transformation is 
an isometry. This assumption is made throughout the article. 

The quantum operation-HTO relation is a generalization of the classical LEP to the 
quantum domain and hence it is primarily of theoretical interest. But, it is possible that 
this may find some applications as well. For example, the restrictions on HTOs can be used 
as a way to theoretically quantify the "information erasing nature" or some other aspect 
of generic quantum operations. In addition to this, the current work can also be useful in 
the design of non-unitary operations (e.g., erasures) on technological implementations of 
quantum devices. If the heat emitted during the operation is desired to be minimized, it is 
necessary to select or adjust the bath levels and the device-bath interaction in a suitable way. 
The explicit construction of the erasure operations that appear in this article is appropriate 
for that optimization task and hence it may be a valuable guide for this purpose. 

The organization of the article is as follows. In Sec. [Tll the HTO associated with a 
realization of a given quantum operation is defined. The nature of the HTO and some of its 
fundamental properties are also discussed in this section. Sec. lIIIl is devoted to the derivation 
of various restrictions on the HTOs. The first one of these is the entropic bound of the 
generalized LEP, which can be expressed for an arbitrary quantum operation. This bound 
is a good starting point for the investigation of the coding-independent characterization of 
quantum operation-HTO relation. After that, quantum operations that completely erase the 
information in the device are taken up and the associated HTOs are described. These results 
are then used for the description of HTOs associated with extremal quantum operations. In 
Sec. llVt some examples are given for quantum operations that do not completely erase the 
information content of the device and it is shown that LEP is not sufficient for a complete 
characterization of HTOs associated with these operations. This indicates the usefulness of 
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the current approach. Finally, Sec. |V] contains brief conclusions. 



II. HEAT TRANSFER OPERATOR 

Let A be a quantum system with a finite number of levels which is used as a memory 
device for holding quantum information. A quantum operation on this system is a trace- 
preserving completely-positive map £ that takes the initial state pa of the device to a final 
state S{pa)- Although all quantum operations are within the scope of this article, special 
attention will be given to those that completely erase all information in the device. The 
operation S will be called a complete erasure if the final state is a constant mixed state po 
for all initial states, i.e.. 



It will be called a Landauer erasure if it is a complete erasure with a pure final state, i.e., 
Po = li'o) (^ol for some \ipo). 

Any quantum operation £^ on A can be realized as a unitary transformation or an isometry 
Uab on a larger system AB, where B denotes an additional system that will henceforth be 
called as the bath. It is assumed that the device and the bath are initially uncorrelated. 
Hence, if the initial state of the bath is pB, the initial state of the composite system is the 
product state pab = Pa® Pb- The final states after the application of the isometry will 
usually be denoted by primes, e.g., p'^^ = Uab{pa ® Pb)U\q. We will say that the bath B, 
the state ps and the isometry Uab form a realization of the operation £ if 



is satisfied for all pA, where here tr^ denotes the partial trace over B. 

It is assumed that the transformation Uab is obtained from the evolution of the composite 
system AB with a time-dependent Hamiltonian HAB{t)- The basic assumption about the 
independence of the process (i.e., the Hamiltonian and therefore the transformation Uab) 
from the information stored in the device is also made. In particular this means that no 
measurements are taken on the composite system AB during the process. Hence, during the 
transformation, any information exchange with A occurs entirely between A and B only and 
this process is described with the Hamiltonian HAB{t)- No third system is involved in the 
transformation except in providing the energy for the necessary work done on AB. 



S{p) = Po for all p 



(3) 




(4) 
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For discussing the energetics of the aforementioned transformation, the time-dependent 
Hamiltonian HAsit) of the composite system AB has to be discussed. The transformation 
process is assumed to be carried out over a finite time interval, say between times tj and 
tf, during which subsystems A and B couple and HAsit) is time dependent. Outside this 
time interval, however, it is assumed that the two subsystems are uncoupled and have time- 
independent Hamiltonians. Thus, we have 



where Ha and if^ are the Hamiltonians of the device before and after the transformation 
respectively and Hb is the bath Hamiltonian. The Hamiltonian of the device is not relevant 
to the main subject of the current article and hence there are no restrictions on Ha and 
H'j^. However, initial and final Hamiltonians of the bath are taken to be identical. This 
is also operationally reasonable for realistic quantum operations where the environment 
corresponds to B; a quantum operation on a device may perhaps structurally change the 
device (so that initial and final Hamiltonians may differ) but it merely couples with the 
environment without changing it structurally. Finally, no assumptions are made about the 
strength of the coupling and the Hamiltonian HAsit) during the time interval ti < t < tf. 

The transformation Uab is the time-development operator corresponding to the Hamilto- 
nian HAB{t) between two times ti and t2 where ti < ti and t/ < ^2- As indicated in Eqs. ([5]) 
and ([6]), it is assumed that A and B do not interact with each other outside the process 
time interval, {ti,tf). For this reason, no special values are needed to be assumed about 
ti^2- However, the results of the current article can also be extended to cases where there 
is a weak coupling between A and B outside the process time interval. In this case, t2 — tf 
should be greater than the thermal-equilibration time scales of the bath. 

Note that for all realistic processes the operator Uab is unitary, i.e., it satisfies U^^Uab = 
UabU\b — ^AB- However, for the derivation of the results in this article, this requirement 
appears to be unnecessarily too restrictive; we only need to assume that Uab is an isometry, 
i.e., only U\^Uab = ^ab is assumed to be satisfied. Therefore, we allow the possibility that 
the transformation Uab is not surjective. The generalization of unitary transformations to 
isometrics also serves another purpose. There are some quantum operations, like Landauer 
erasures, which can only be realized by proper isometrics instead of unitaries. This is quite 



HABit) = H'^ + Hb ioit>tf , 



HABit) = Ha + Hb for t <ti and 
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obvious for Landauer erasures since all possible accessible final states do not span the whole 
of the Hilbert space of AB, because any initial pure state of AB is mapped to a final product 
state where A is in \ipo) state[24']. As Uab is not surjective, it cannot possibly be unitary. 
Such proper isometries can in principle be realized by a Hamiltonian evolution, provided 
that ideal tools like perfectly impenetrable barriers are used for making parts of the Hilbert 
space inaccessible. 

It is assumed that the bath is in canonical thermal equilibrium state before the transfor- 
mation and thus 

Pb = ^exp (-/^iJs) , (7) 

where /3 = l/ksT is the inverse temperature and Zb is the corresponding partition function. 
It is appropriate to also assume that the bath is sufficiently small so that is manageable. 
We will say that B is a finite system if Zb is finite for all temperatures {Zb 7^ 00). In 
this case, B has finite thermodynamical properties. Note that a single Hydrogen atom or a 
single free particle placed in an infinite space has an infinite partition function and therefore 
these are not finite systems by the current definition. But, if the atom or the particle is 
placed inside a large but finite box, then they will become finite systems. Note also that any 
realistic system, which has a finite number of particles constrained inside a finite volume, 
will be classified as finite by this definition. Infinite systems have individually vanishing 
level occupation probabilities and consequently the analysis of such systems present some 
difficulties. This is the main reason for using finite systems in this article. 

Let Eb = ti pbHb be the equilibrium energy of the bath. Since there is no structural 
change in the bath, the heat dumped to the bath during the realization is equal to the 
change in the average energy, AEb, between times ti and ^2- It can be expressed as 

AEb = ti{p'^ - pb)Hb , (8) 
= tiAB (pA (S) Pb) (u^Bi^A ® Hb)Uab^ " Eb 
= tr^ (paQ) 

where Q is the HTO of this realization, which is defined as 

Q = tiB [{1a ® Pb) U\b{1a ® HbWab] - EbIa • (9) 

A number of comments have to be made in order to comprehend the real nature of the 
HTO. First, Q is a hermitian operator defined on the state space T-La of the device A. Its 



expectation value is evaluated with the initial state, pA, of the device. In other words, it is 
an operator that can be used for expressing the energy transferred to B in terms of the state 
of A. Second, only the average of Q can be given a physical meaning. It has been shown 
that no quantum observable can be defined for the work done or the heat transferred during 
a process that captures the detailed distribution of these quantities |2Q.]. However, if only 
the average Ai?^ is needed, then Eq. (|8]) is the correct definition, which eventually yields 
the HTO Q. On the other hand, Q cannot be used for computing the higher moments of 
the energy transferred to B. For example, even though the average energy transferred to B 
is given by tr paQ, the quantity tr paQ^ or the averages of other nonlinear functions of Q do 
not have the anticipated meaning. Averages of such nonlinear functions of energy change of 
B should be computed separately following the approach in Ref. |20i. This crucial property 
of the HTO must always be kept in mind. 

Once the HTO is computed, the average total work that needs to be done for carrying 
out the operation can be obtained as 

W^tot = tr {£{pa)H'^ - PaHa + PaQ) • (10) 

As a result, for investigating the average energies transferred during the realization of the 
operation S, the only nontrivial quantity to be determined is the heat dumped to the bath. 
This is the main reason why this article concentrates solely on the HTO. 

The isometry condition U^qUab = ^ab and the fact that the bath is in thermal equilib- 
rium imposes some restrictions on the operator Q. The investigation of these restrictions 
is the primary subject of this article. The main question that needs to be answered is the 
following: for any given quantum operation hermitian operator Q and temperature value 
T, is it possible to find a bath B (namely a Hamiltonian Hb) and an isometry Uab such 
that both Eq. (jlj) and Eq. ([9]) are satisfied? We will say that Q is a possible HTO for £ 
when such a bath and isometry exist. Our main problem is then to find all necessary and 
sufficient conditions for Q to be a possible HTO for £ at temperature T. These conditions 
will enable us to discuss the thermodynamics of quantum information processing (for given 
£) without specifying anything about the bath, the bath-device coupling and the specific 
way the quantum operation is realized. 

Before embarking on the investigation of this problem, it will be convenient to list the 
following basic properties of the HTOs. 
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(i) First, note that any restriction that the HTOs obey can be expressed entirely in terms 
of £ and the ratio Q/ksT = /3Q. To see this, note that both S and PQ can be expressed 
exclusively in terms of Uab and pb, with no explicit T (or Hb) dependence, i.e., Eq. (jlj) and 



f3Q 



tlB {1a® pB)Uls{lA^lnpB)UAB 



+ lAtTB{PB Inps) 



(11) 



Hence, the main problem becomes: for given 8 and /3Q, can we find an appropriate pb and 
Uab such that Eqs. (jlj) and ffTTl) are satisfied? 

(ii) If Q is a possible HTO for an operation S, then for any positive number a, the 
operator Q = Q + qIa is also a possible HTO for the same operation. In other words, it is 
always possible to dump an additional heat a whose amount is independent of the state of 
the device. 

To show this, suppose that Uab on bath B is a realization of S that yields the HTO Q. Let 
C be a two-level system with Hamiltonian He = Aaz/2 where denotes the associated Pauli 
spin matrix. Now, take BC as the new bath and apply the isometry Uabc = Uab ® {o'x)c to 
the composite system ABC. Since the transformation that C is subjected to is independent 
of how AB is changed, it is obvious that Uabc realizes the same operation S. But, the heat 
dumped to the bath BC is now 



where AEq = A ta,nh.{A/ 2k bT). If the value of A is chosen such that AEc = a, then we 
get AEbc = ^"i^PaQ- This shows that Q is also a possible HTO for the same operation. 

(iii) The set of two-tuples {£, Q) formed from quantum operations £ and allowed HTOs Q 
is a convex set. In other words, if Qi is a possible HTO for the operation £i for i = 1, 2, . . . , n, 
then for any Aj > with Aj = 1, the operator Q = J2i ^iQi is a possible HTO for the 
operation £ = J2i ^i^i- 

Before showing this, first note that realizations for each £i can be assumed to involve the 
same bath B. This can always be done by taking the bath sufficiently large. As a result, 
without loss of generality, suppose that U% is an isometry on AB that realizes £i and yields 
the associated HTO Qi. Let C be an n-level quantum system and let {|1) , ■ ■ ■ , |n)} be an 
orthonormal basis in the Hilbert space of C. Let the initial state of C be pc = Xli K) (^1- 



AE^ 



BC — 



tiipAQ) + AE, 



•c 



(12) 
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Now, consider the isometry Vabc on ABC given by 

n 

VABC = Y.UfB^i\^){^\)c (13) 
i=l 

It is a simple exercise to show that this isometry gives rise to the desired convex combinations, 
i.e., Vabc reahzes the operation £ = J2i '^i^i the associated HTO is Q = X]j '^iQi- Note 
that, in here, C acts hke a controller system for deciding which operation should be applied 
on AB. The decision is given probabilistically such that £i is applied with probability A,. 
Note also that the average energy of C does not change during this process. In such a case, 
the overall operation and the associated HTO are given by the expected convex combination 
expressions. 

(iv) As a special case of the convexity result in (iii), the set of allowed HTOs associated 
with a fixed quantum operation £ is also convex. This set is not bounded from above by 
(ii); but it is bounded from below. We are primarily interested in finding lower bounds, if 
possible tight ones, for this set. 



III. RESTRICTIONS ON HEAT TRANSFER OPERATORS 



A. Entropic Restriction of LEP 

Before investigating the relations between HTOs and quantum operations, it is appro- 
priate to start with the state-dependent entropic restriction of LEP. Let S{p) = — trplnp 
denote the von Neumann entropy of the density matrix p. For an arbitrary quantum oper- 
ation, the average heat emitted to the bath is bounded from below by the drop in the von 
Neumann entropy of the state of the device. 

Theorem 1 (Generalized LEP) If Q is a possible HTO for an operation £, then 

tl PAQ>kBT[S{pA)-S{£{pA))] , (14) 

for all states pA of the device A. For a given pa, the equality sign in holds if and only 
if the relation p'j^^ = £{pa) ® Pb is satisfied for some (and therefore for all) realization that 
yields £ and Q. 

Note that the inequality above essentially gives infinitely many bounds on the HTO, 
i.e., there is a separate bound on Q for each possible device state pa- The generalized 
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14j . especially in the context of classical 



LEP has been discussed in several works 
information processing. However, for the purpose of treating the equality case, the following 
straightforward derivation is included in here. It can be deduced simply from the following 
wonderful identity for deviations from canonical thermal equilibrium states 

AEb = k^T {S{p'b\\pb) + S{p'b) - S{pb)) (15) 

where 

^((t||p) = tra(lna-lnp) (16) 

denotes the relative entropy function. Only two inequalities will be used for converting 
this identity to the relation in Eq. f|T^ . The first one is the nonnegativity of the relative 
entropy: S{(r\\p) > with equality holding if and only if a = p. The second inequality is the 
subadditivity of the von Neumann entropy for the final state, i.e., ^'(p^^) < S{p'^) + S{p'b) 
where the equality holds if and only if the state is in product form p^^ = p'^® p'^- Finally, 
using the fact that the initial and final states are related by an isometry we conclude that 
p'j^B ^'^d pab are isospectral. In particular, S{p'j^b) = S{pab) = S{pa) + S{pb) and hence 

AEb > kBT{S{p'B) - S{pb)) > kBT{S{pA) - S{p'^)) , 

which is the desired inequality. It is easy to see that both of the inequahties become equalities 
if and only if p^^ = S{pa) ® ps.D 

Note that the equality can be satisfied only in some exceptional cases. When that hap- 
pens, we have p^^ = £{pa) ® Pb and therefore the state of the bath is unchanged, i.e., 
p'b = Pb- Consequently, the average heat dumped to bath is ti paQ = 0. Moreover, since 
the initial and final density matrices of AB are isospectral, it follows that pA and £{pa) must 
also be isospectral. 

It is of some interest to investigate the issue of whether the bound in ( JT^ can be improved 
by considering the device A to be part of a larger system AX where X is not involved in the 
realization. In other words, the inequality 

tr p^Q > kBT{S{pAx) - S{{£®X){pAx))) (17) 

is also valid for any X and any initial state pax- Can this generalized inequality be tighter 
than Eq. (fT4l) ? The answer is negative. To show this, introduce another system Y such 
that AXY is a purification of Pax- The isometry of the realization is applied on AB and 
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hence the system XY is untouched. Using primes for denoting the final states we have 
S{pa) = S{pxy) = S{p'xy) and S{pax) = S{py) = 5'(py). Then we invoke the following 
alternative form of the strong subadditivity of the von Neumann entropy 21 1 



Sip'xv) + S{p'xa) > Sip'y) + Sip'^) 



(18) 



to arrive at 



S{pa) - S{p^) > S{pAx) - S{pax) 



(19) 



The last inequality clearly shows that the lower bound given by Eq. f[T7l) is always smaller 
than the bound in f|T^ . Hence, the inclusion of additional systems cannot possibly make an 
improvement in the entropic bounds of LEP. 

B. Restrictions for Complete Erasure Operations 

Even though Theorem [1] places some restrictions on the possible HTOs, it is possible 
that these may not be complete. In other words, there might be operators Q which satisfy 
the inequality f[T^ for all possible initial states pA, but it still may not be a possible HTO. 
Some examples of such operators will be discussed in Sec. El However, it appears that 
the inequalities (HM are sufficient when £^ is a complete erasure operation. This subsection 
contains the proof of this statement. 

But, before that, a concise form of the inequalities in f|T^ should be found and for 
this reason it is necessary to introduce a real-valued function of hermitian operators. The 
needed function is the Legendre transform J of the von Neumann entropy function S. For 
any hermitian G, this function is defined by 



In order to arrive at the main property of this function, first note that for any hermitian 
operator G, the operator 



J{G) = - In (tre-^) 



(20) 



ac = exp(J(G'))e 



-G 



(21) 



is a density matrix. Then, using the nonnegativity of S{p\\aG) we find 



J{G) + S{p) <tipG 



(22) 
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an inequality which is vahd for any density matrix p and any hermitian G. From this 
inequahty it is possible to derive the following Legendre transformation identities 

J(G) = min(trpG-5(p)) , (23) 

p 

S{p) = inf (tr pG-J(G)) , (24) 

where in Eq. ( 123|) the minimization is over all density matrices and the minimum is reached 
when p = (Jg- In Eq. f l2^ the infimum is over all hermitian operators G. If p has no zero 
eigenvalues, then right-hand side of that equation is minimized by G = — In p. Otherwise, 
if p has zero eigenvalues, then the right-hand side of f l2^ has no minimum; but infimum is 
obtained by a sequence of density matrices which approach to the limit = — In p, which 
is now an operator that has some infinite eigenvalues. 

Consider a complete erasure operation £ where £{pa) = Po for all states pA- In this case, 
the inequality (fT^ . which is valid for all density matrices pA, can be written compactly 
using the J-function as 

Jim>-S{po) . (25) 

Moreover, this condition is almost sufficient for Q to be a possible HTO. There is only a 
slight complication in the treatment of the equality case which is related to the finiteness 
property of the bath. The following theorem gives the complete characterization of HTOs 
for complete erasures. 

Theorem 2 Let Q be a hermitian operator and £ be a complete erasure to the final state 
Po . When only realizations involving finite baths are considered, Q is a possible HTO for £ 
if and only if 

(z) either J{l3Q) > -S{po), 

(a) or J{f3Q) = —S{po) and a^q is isospectral with pg. 

Therefore the equality sign in fl25l) applies only under very rare situations where the 
eigenvalue spectrum of Q must satisfy strict conditions; in particular the dimension of the 
state space of A must be equal to the matrix rank of the final state po. If this is not the 
case, then the strict inequality holds. 

The proof of the necessity of the conditions (i) and (ii) has already been argued above 
in obtaining Eq. (|25|1 : only a treatment of the equality case is left. If equality holds in 
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(125|) . then the density matrix that minimizes the expression (123|) . namely pA = (^PQ-, gives 
equahty in f|T4|) . and therefore by Theorem [H pA must be isospectral with £{pa) = po, i.e., 
the condition in (ii) holds. This completes the proof of the necessity of the conditions of 
Theorem [21 

For the sufficiency proof, we have to start with a given Q that satisfies either (i) and 
(ii) and then show that a realization that yields S and Q can be constructed. The needed 
isometry is actually quite simple. It essentially copies the entire information on the device 
to some subsystem of the bath, so that the information is completely erased, and at the 
same time it copies the thermal equilibrium state of some subsystem of B to the device, so 
that the device has always the same final state. Even though the isometry is simple, the job 
of adjusting the HTO to be equal to the given operator Q is technically complicated. For 
this reason, the proof of the sufficiency of the conditions in Theorem |2] is done separately 
in Appendix [Bl However, it is appropriate to give a separate proof in here for the special 
case of Landauer erasures for which the required bath is simpler. Treatment of the technical 
details in this proof will also simplify the presentation of the full proof of Theorem [2] in 
Appendix [Bl 

Let £^ be a Landauer erasure to the standard state \ipo) {£{pa) = iV'o) (V'ol for all Pa)- Note 
that the special case (ii) of Theorem [2] can hold only if the dimension of the Hilbert space 
of A is (i = 1. This is a very trivial case where the system A cannot hold any information 
at all. In such a case, the HTO Q is a number and the theorem just says that this number 
is nonnegative, Q > 0. For the generic case of (i > 2, however, only the strict inequality 
of (i) holds. The following is merely a corollary of Theorem |2] but we state it as a separate 
theorem and give a separate proof. 

Theorem 3 Let £ he a Landauer erasure on a device A having d levels with d > 2 and let 
Q be a hermitian operator. The following are equivalent. 

(i) Q is a possible HTO for a realization on a finite bath. 

(ii) tr e~^*^ < 1. 

(Hi) ti paQ > kBTS{pA) for all pA- 

It is obvious that the statements (ii) and (iii) are equivalent; both are simply stating 
that J{fiQ) is a strictly positive number. The statement (iii) is the conventional way of 
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Bath 



FIG. 1: The pictorial description of the isometry used in the proof of Theorem [2] for a Landauer 
erasure operation. Here, the state of the device A is made equal to the standard state {ipo) while 
the original information content of each subsystem is shifted to the right. 

expressing the LEP. Statement (ii) is the quantum analog of the inequahty in Eq. ([1]). (The 
strict inequahty is a minor mathematical detail; it arises from our insistence on using finite 
baths.) 

It has been shown in Theorem [1] that (i) implies (ii) and (iii). What is left is the proof of 
the implication (ii)^(i). Let Q be a hermitian operator such that J{PQ) > 0. A realization 
of the Landauer erasure will be constructed below and it will be shown that the associated 
HTO is equal to Q. 

a. The bath and its basis states. It is obvious that the bath needed for this purpose has 
infinitely many levels. Only in that case, it is possible to copy the information on A onto 
B and at the same time, retain the information already present in B. For better visualizing 
the information flow, we will think of the bath as being composed of infinitely many copies 
of a d-leve\ system X. The constitution of B can be expressed as B = X1X2X3 ■ ■ ■ . Let the 
standard orthonormal basis of X be denoted by {\i)} where i = 0,1, . . . ,d — l. The following 
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self-explanatory shorthand will be used for the state vectors of the bath 

I«i,«2,i3, = ® N2)x2 N3)x3 ■ (26) 

The Hilbert space "Hb of the bath will be defined as the linear span of the above states where 
the sequence {ik} contains only finitely m,any nonzero elements. When the Hamiltonian is 
defined below, it will be seen that the states which are not included into Hb have infinite 
energy and zero equilibrium-state probabilities. Consequently such states will never enter 
into the discussion. 

In addition to this, defining T-Lb in such a way makes it a separable Hilbert space having 

a countable orthonormal basis. To see this, note that the sequence {ik} can be considered as 
the digits of the base-d representation of an integer n, written in short as n = (• • ■izi2ii)d- 
The integer n is given in terms of the sequence as 

oo 

n = ii + i2d + i^d? ^ ^^ikd!"'^ . (27) 

fc=i 

Consequently, there is a one-to-one correspondence between nonnegative integers n and 
sequences {ik} that has finitely many nonzero elements. For this reason, the standard 
orthonormal basis oH-Lb can be enumerated by a nonnegative integer n and therefore the 
alternative label 

= Ni,«2,«3,---)b (28) 

can be used. It is thus clear that the Hilbert space 1-Lb has a countable basis and hence it 
is separable. 

h. The isometry. Let Qi {i — 0,1, . . . ,d—l) he the eigenvalues of Q arranged in nonde- 
creasing order and let {\i)} be the orthonormal basis of T-La formed from the corresponding 
eigenvectors. 

d-l 

<9 = I]9^N)(^I {qo<qi<---<qd-i)- (29) 

1=0 

The isometry Uab on the composite system AB is defined as follows 

C/ab|^)a® Ni>^2,«3,---)b = |V'o)a® K>^i>^2,---)b , (30) 

for any allowed values of input labels. In words, the isometry sets the state of A to the 
standard state j'^o) and copies its contents onto the subsystem Xi. It also forms an infinite 
chain of information fiow inside the bath as the contents of Xj. are copied onto X^+i {k > 1) 
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for all k. The flow of information is pictorially depicted in Fig. [TJ Essentially the isometry 
shifts the information content of all subsystems of AB onto the next subsystem on the right. 

This flow is highly reminiscent of the motion of the guests in Hilbert's Hotel Infinity 
paradox. In that hotel, all of the infinitely many rooms are occupied. But, it is still possible 
to accommodate a new guest by moving everyone to the next room. The isometry Uab 
essentially does the same for information. It retains the whole information already stored 
in B while it squeezes additional information coming from A into B. Our hotel, however, 
has an important difference from Hilbert's hotel: as it follows from LEP, there is a nonzero 
check- in cost. In other words, while squeezing new information into B, the average energy of 
B increases. To be able to compute that increase, the bath Hamiltonian must be specified. 

c. The Hamiltonian of the bath. The Hamiltonian will be constructed such that each 
subsystem of the bath is independent and the energy eigenstates are identical with the 
standard basis states. The precise energy eigenvalues of the subsystem Xi will be chosen to 
be identical with those of Q (plus a possible shift), i.e., Hx^ = Q + cilxi in the standard 
basis where here, ci is some constant. The energy eigenvalues of the remaining subsystems 
X2, X3, ... will only shift the HTO by a multiple of identity. With this choice, the energy 
change of the bath is AEb = tr(p^(5) + C2 where C2 is a constant that is independent of 
Pa- The value of C2 depends on the energy levels of X2, X3, . . .. By adjusting these levels 
appropriately, it is possible set C2 = 0. This shows that the HTO is Q and completes 
the proof. Since the details are complicated, the adjustment of the energy levels of other 
subsystems is discussed in Appendix |Al □ 

The realization described in the proof above can be used for finding conditions for min- 
imizing the average heat emission during a Landauer erasure operation. Suppose that A 
is used as a memory device in an implementation of a quantum computer. Suppose also 
that A is subjected to a Landauer erasure at ambient temperature Tenv in some steps of 
the operation of the computer. Let denote the average state of the device before the 
erasures are applied. The average heat emitted per erasure is then g = tr p^Q where Q is 
the HTO. The problem is to minimize q under the restriction that Q satisfies the conditions 
mentioned in Theorem |3l According to Theorem [H q is bounded below by kBTcmS{p^^). A 
simple inspection shows us that this bound is attained by a small error by the HTO 

Q = -^BTenvlnpX + elA (31) 
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where e is a small positive number. The construction in the proof of Theorem [3] can then 
be followed for adjusting the energy levels of the bath, the detailed Hamiltonian HAB{t) 
during device-bath interaction and the consequent isometry. In a realistic situation where 
the environment serves as the bath B, it is not possible to modify the energy levels of B at 
will and the Landauer bound may not be attained; but the construction above can still be 
useful. For example, a finite chain X1X2 ■ ■ -Xtv can be used for connecting the device and 
the environment. By adjusting the energy levels of the subsystems in the chain and using 
the same information flow, q can still be significantly decreased. 



C. General Quantum Operations 

A generic quantum operation causes only a partial erasure of the quantum information 
stored by the device by forming an entanglement between the device and the bath during 
the realization. Even though the inequality ( |T^ continues to be valid for the generic case, 
it does not necessarily give a complete characterization of the possible HTOs as it will be 
discussed in Sec. |Vl This subsection contains a few basic results about the HTOs of generic 
quantum operations that surpass (|T^ in implications. 

Let £^ be a quantum operation on the device A. Such an operation can be written in the 
Kraus form 

n 

8{pA) = Y,^ipM (32) 

i=l 

where the Kraus operators Mj satisfy 

n 

Y,MlMi = lA . (33) 

i=l 



Kraus operators can be chosen in various different ways |2 II]. However, a representation that 
has the smallest number n of terms will be preferred in this article. In this case, the set of 
Kraus operators {Mi, . . . , M„} are linearly independent. 

Now, let Uab be the isometry for a realization of the operation £. This isometry can be 
written as 

n 

Uab = J2^Mi)A ® {L,)b (34) 

i=l 

for some operators Lj acting on the Hilbert space of the bath B. This can be shown by 
completing the set of operators {Mi, M2, . . . , M^} to a basis for the operator space by 
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adding new, linearly-independent operators M„+i,M„+2, ■ ■ ■ to the list, and then applying 
the expansion in Eq. f l34|) with a larger upper limit of summation. As Uab is a realization 
of we have 

£{Pa) = tlBUAB^pA® Pb)U\b (35) 
= M.,paM] (tr PbL]L^ . (36) 

ij 

The linear independence of {Mi} implies the linear independence of the set of maps 
{MiPAMj}ij on the density matrices. For this reason, we have 

, 6ii for i,j < n , ^ 

tTpBL]L,= \ (37) 
I for i > n OT j > n 

and consequently ti pBL\Li = for all i > n. Since pb is a canonical equilibrium state, it 
has full rank and therefore = for all i > n. This shows that the summation in Eq. (IMll 
contains only n terms. The following theorem directly follows from this expansion. 

Theorem 4 Let E he the quantum operation given in / f3^) where Mj are linearly indepen- 
dent. 

(a) If Q is an HTO for £ , then there is an n ^ n hermitian matrix q, (which might he 
called as the heat transfer matrix^ such that 

n 

Q=Y,%MlM, . (38) 

(h) If n > 2 and the set ofn"^ operators {mJ Mj}^j^-^ is linearly independent, then Q is an 
HTO for S if and only if the unique matrix q in / fS'Sj) satisfies 

tre-«/'=«^<l . (39) 

For n = 1, Q is an HTO if and only if Q = at a where a > 0. 

(c) If {mJ Mj}'^j^-^ is linearly dependent, then the condition is sufficient for Q given 
hy (E^j to he a possihle HTO; hut that condition is not necessary. 



Note that the set of possible quantum operations on A is a convex set. It can be shown 
that the extreme points of this convex set correspond to quantum operations for which the 
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set of operators {M^ Mj}^j^^ is linearly independent 22 1. Therefore, part (b) of the theorem 
above gives a complete characterization of the HTOs of all extreme quantum operations. 

Proof: For (a), the use of the definition of the HTO given in Eq. (jH]) directly leads to 
Eq. (138 p where 



Qij = tr {^pbL\HbL^ - Es^ij . (40) 
To show (b), first write the isometry condition U\qUab = Iab as 

n 

Y,{MlM,)A<^{L\L^)B = tA<^tB . (41) 

The linear independence of {mJ Mj}^j^i then leads to the following relations 

LlL,=6ijlB . (42) 

In other words, each Lj is an isometry on the bath that maps Hb onto mutually perpendicular 
subspaces of Hb- The last set of relations is actually equivalent to the statement that Uab 
is an isometry. 

Consider now a hypothetical n-level system R. Let be an orthonormal basis for 

the Hilbert space of R. The following operator on the composite system RB, 

n 

W^i?B = $^(|l)(^|)«®(^.)B (43) 

i=l 

is then an isometry. It can also be readily verified that Wrb is an isometry if and only if 
Uab is also an isometry. The most important point in here is that Wrb is a realization of a 
Landauer erasure on R. Furthermore, the HTO associated with this erasure is given by 

Q'^ = trs [(1r ® pb) W^si'^R ® Hb)Wrb] - EbIr . (44) 

id 

For n > 2, the application of Theorem [3] proves the claim. The result for the special case of 
n = 1 also follows trivially from the discussion in this subsection. 

Part (c) follows from the converse of the argument used above for part (b). Suppose 
that q satisfies (l39ll . In that case, Q'p, given by Eq. (115|) is an allowed HTO for a Landauer 
erasure. Hence one can construct the associated realization; let Wrb be the isometry for 
this realization. We use Eq. fH3|) to define Lj and then use Eq. fl34l) to define Uab- It is 
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straightforward to verify that Uab is the desired isometry that reahzes S and has q as the 
heat transfer matrix. 

Finally note that any q that satisfies fl39p necessarily yields a positive definite Q. However, 
Theorem [2] allows HTOs that may have negative eigenvalues. This shows that the condition 
(I39p is not necessary. The proof of Theorem H] is thus completed. □ 

Note that part (b) of Theorem H] is proved by finding a closely related Landauer erasure 
on a hypothetical system. It is possible to extend the scope of part (b) by making use of 
complete erasures. However, the expressions used in such an extension do not appear to be 
particularly illuminating and hence they are not given in here. Theorem S] is a good starting 
point for discussing the properties of HTOs associated with generic quantum operations. 
The following result is an example. It is essentially a different way of saying that the set of 
HTOs associated with a given quantum operation has no upper limit. 

Theorem 5 Let q be a heat transfer matrix for £ which is related to the corresponding HTO 
by Then, for any q' with q' > q, q' is also a possible heat transfer matrix. 

To prove this, first define the matrix s = q' — q and observe that it is strictly positive 
definite, i.e., all eigenvalues of s are positive. In such a case, it is possible to find a positive 
number t > 1 such that the matrix q + ts satisfies ( l39l) and therefore g + ts is a possible heat 
transfer matrix by part (c) of Theorem |H Finally, by the convexity of the set of HTOs, the 
matrix 

q' = ]{q + ts) + ^-^q (46) 

is a possible heat transfer matrix. This completes the proof. □ 

Note that this theorem cannot be extended by replacing the matrices q by the associated 
HTOs mainly due to the restriction in Theorem H] (a) which must be satisfied by all HTOs. 
It is plausible that the strict inequality q' > q can be replaced by the weaker requirement of 
q' > q, but this conjecture is still awaiting proof. 

IV. DIFFERENCES FROM LANDAUER'S BOUND 

The restrictions satisfied by the HTOs are closely related to LEP since both follow only 
from the condition that the transformation on the composite system is an isometry. However, 
LEP only relates the dumped heat to the drop in the von Neumann entropy of the state of the 
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device. For complete erasures, it turned out that LEP also gives a complete characterization 
of the HTOs. In other words, Theorems [2] and [3] are identical in content with Theorem [H 
they are just expressing the same statement using different equations. As it will be shown 
below, this is not the case for other quantum operations, i.e., it is possible to find examples 
of operations £ that are not complete erasures, such that the HTOs associated with £ 
satisfy further restrictions that cannot be explained by LEP alone. Such operations can be 
considered to be "partial erasures" as some part of the initial information present in the 
device will be retained (since £{p) depends on p). The Landauer bound in Theorem [1] still 
applies to such £, but this bound fails to imply the additional features of HTOs that are 
stated in Theorem |H 

For example, the special structure of the HTOs given in Theorem H] (a) does not appear 
in the generalized LEP inequality f|T^ . To see this clearly, consider unitary rotations (or 
isometrics) as quantum operations, i.e., £{pa) = W paW^ where W'^W = Ia- In such a case, 
the entropic bound f fT4l) merely states that the HTO must be a positive-semidefinite operator. 
However, in reality, the HTOs should not only be positive-semidefinite but they should also 
be a multiple of identity as asserted in Theorem H] (b). Only in such a case, the average 
heat dumped to the bath becomes independent of the state of the device. Otherwise, any 
Q which is not a multiple of identity would create serious problems for quantum mechanics: 
it implies that by investigating the energy of the bath, some information can be obtained 
about the state of A without disturbing it! It is apparent that Theorem H] (a) is intimately 
related with the characteristic features of quantum information. 

Apart from this, there are situations where a given Q satisfies the special form given in 
Theorem H] (a) and the entropic inequalities ffT^ . but it still fails to be an HTO. As an 
example, consider the special problem of finding an HTO of the form Q = oIa in such a 
way that a has the minimum possible value. If the quantum operation is extremal, then by 
Theorem m (b), Q is an HTO if and only if a > kBTlnn where n is the minimum number 
of Kraus operators (for n > 2). However, if only the entropic inequalities (HM are used, it is 
possible to obtain a different lower bound for a, especially if the quantum operation is near 
to identity. 

As an example, consider the one-parameter family £t of quantum operations defined by 
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n = 2 Kraus operators 



Mi(t) = tX , (47) 
M2(t) = - t^X^X , (48) 

2 

^t{pA) = 5^M,(t)tp^M,(t) , (49) 

i=l 

where X is a non-hermitian operator such that {M^ {t)Mj(t)}'^ are hnearly independent. 
Note that in the hmit t — y 0, the operation St approaches to the identity transformation. 
Therefore, the maximum drop in the von Neumann entropy of the device, 

Bt = m^x{S{pA)-S{St{pA))) , (50) 

PA 

approaches zero in the same hmit. Then, the smaUest possible value of the coefficient a 
allowed by the entropic restrictions of LEP is amin = ksTBt, which also tends to zero as 
t — > 0. However, it is shown above that a > ksT ln2, independent of how close St to the 
identity operation. This clearly shows that Theorem H] (b) contains restrictions on HTOs 
that are not implied by the entropic LEP restrictions. 

In order to eliminate a possible misunderstanding that may arise from the previous ex- 
ample, let us stress on the fact that once we lift the restriction that Q is proportional to the 
identity, it is possible to find HTOs that approach to the zero operator when St approaches 
to the identity operation. For example, for St given above, 

2 

Q = J2<liMi{tyM,{t) , (51) 

i=l 

with (3q2 = and ^ ln(l/t2) yields an HTO Q with 

PQ ^ [Ia + Hl/e)X^X) (52) 
which converges to zero as t — )■ 0. 

V. CONCLUSIONS 

In this article, the concept of heat transfer operator is introduced and some of its prop- 
erties are discussed. The HTO is the main quantity needed for analyzing the energetics of 
quantum operations. Hence, the complete characterization of the HTOs associated with a 
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given operation 8 occupies a central place in the thermodynamics of quantum information 
processing. This article mainly contains results on this problem. A complete characterization 
could be provided only for complete erasures and extreme operations. The characterization 
of HTOs of generic non-extremal quantum operations is still an open problem. 

The set of HTOs associated with a given quantum operation might be considered as 
a means of describing the disturbance caused by the operation on the environment. This 
description is provided through a set of operators that has a clear operational interpretation. 
This points out the main theoretical interest for this problem. In addition to this, both the 
results obtained and the techniques used in this article can be useful in the implementation 
of non-unitary operations in future quantum computers. For example, the minimization of 
the average heat emission during an erasure operation requires a careful construction of the 
transformation on the extended system of the device and the bath. The constructions used 
in this article will be useful for such design problems. 
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Appendix A: Adjusting Bath Levels in the Proof of Theorem [3] 

It is convenient to think of the subsystem Hamiltonians to be dependent on a continuous 
parameter v. For this purpose, consider d continuous functions ej(z/) (i = 0, 1, ■ ■ ■ , (i — 1) of 
a real parameter v G [0, 1], which will be used for energy eigenvalues of subsystems. They 
are chosen to satisfy the following properties: 

(i) The zeroth level has zero energy, eo(z/) = 0, for all v. 

(ii) The energies of all excited levels diverge to +oo as v approaches to 1 (i.e., 
lim,^^i ej(z/) = oo for any i > 1.) 

(iii) The energy spectrum at z/ = is identical with the eigenvalues of Q up to a constant 
shift as ej(0) = qi — qo for all i. 

Let {z/fc} be an infinite, increasing sequence starting from and converging to 1, i.e., 
1^1 = < 1^2 < lys < ■ ■ ■ and limfe_j.oo i^k = 1- The energy levels of will be taken as eiivk)- 
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In other words, in terms of 

d-l 

h{u) = Y,^,{u)\t){t\ , (Al) 

1=0 

the Hamiltonians of individual subsystems can be expressed as Hx^. — h{v}.). The bath 
Hamiltonian Hb is the sum of the subsystem Hamiltonians and therefore the energy eigen- 
value of the bath state \n) ^ = \ii,i2, is, - ' ' ) b 

oo 
k=l 

Note that if {ik} had infinitely many nonzero entries, then the corresponding energy would 
be infinite. This justifies the definition of the Hilbert space "Hb in the way described above. 

Let (7(i/) = exp(— be the parameter-dependent density matrix and (^{v) be 
the corresponding partition function. The shorthand notation 0"^ = (j(^'fc) will be used for 
the thermal equilibrium states of X^. The partition function of the bath is then given by 

oo 

ZB^Jla^k) . (A3) 

k=l 

Note that the convergence of this product depends only on the large k behavior of the 
sequence {fk}- By adjusting the way Vk converges to 1, it is possible obtain a finite Zb- 

Let us now compute the HTO. In the following expressions, the subscripts are used for 
indicating the subsystem a given density matrix applies to. Suppose that the initial state of 
the device A is p. The initial and final states of AB are given by the following self-explanatory 
expressions. 

pAB ^ { P )a® ((X2)X2 (A4) 

p'ab = (iV'o) (V'oI)a ®{p)x,® ® • • • ■ (A5) 

The change in the average energy of the bath is given by 

A£;b = tr(pg) - kBTJ + A (A6) 

where J = J{(3Q) and 

oo 

A = kBTJ - go - tr CTiifxi + ^ tr(cTfc-i - (JkjHxu , (A7) 

k=2 
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which is a constant independent of the device state p. Let us first note that, by using the 
identity in Eq. f fTSj) . it is possible to express A as 

oo 

^ = kBTY,S{ok\W+i) , (A8) 
fc=i 

which clearly shows that A is a strictly positive number. It will be argued below that the 
sequence {uk} can be selected in such a way that A can be adjusted to be equal to any 
positive number and this can be achieved with a finite Zb- This will be done by showing 
that A can be made arbitrarily large and arbitrarily small and then allude to continuity to 
prove the claim. 

(1) First, note that if the parameter V2 of X2 is increased towards 1, the final energy of 
X2 approaches to infinity. Hence, in this limit A diverges to 00. Therefore, A can be made 
arbitrarily large. (2) In order to make A small, the parameters of neighboring subsystems 
should be chosen nearly equal. In other words, the increments Vk+i — i^k should tend to 0. 
In this limit, the series in Eq. (]A7p approach to an integral. 



Jo du 

= ksTSiaiO)) . (A9) 

Using h{0) = Q — QqI, it can be seen that A approaches to in this limit. This shows that 
A can be chosen to be as small as possible. Moreover, for any given parameter sequence 
{z/fc} it is possible to adjust the large k behavior of these sequences and make Zb finite. For 
example, such an adjustment can be made only for indices k that exceed a threshold /cmin- 
By selecting the threshold k^i^ sufficiently large, it is possible to insure that A does not 
change much. In summary, therefore, A can be made as small as possible by keeping Zb 
finite. 

Since A depends continuously on u^, it is then possible to adjust these parameter se- 
quences such that A has the value A = fc^TJ, which is a strictly positive quantity by the 
initial supposition. Then, by Eq. flA6p . the change in the energy of the bath is AEb = tr pQ 
and hence Q is the HTO associated with this realization. This completes the proof of the 
Theorem [31 
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Note that if the finiteness restriction on the bath is relaxed, then the strict inequahties of 
(ii) and (iii) of Theorem [3] will be a non-strict inequality. The construction of the realization 
above can also be adapted to the special case of J{/3Q) = 0. In this case, all subsystems of 
B are in the same ai state and Zb is necessarily infinite. Unfortunately, the correct energy 
change of B cannot be calculated directly in such a state. At first sight, it appears that 
the isometry changes the state of only the subsystem Xi of B. From here, one may try to 
compute the change in energy of B as AEb = tr(p — ai)Hxi, which can easily be seen to 
be incorrect. There is a subtle energy contribution of the shift of information in an infinite 
chain and that contribution can only be computed correctly as in the proof above by using 
a bath with finite Zb- 

Appendix B: Sufficiency Proof of Theorem [2] 

This section contains the proof of the sufficiency of the conditions in Theorem O A 
separate proof has to be given for the two conditions. For the special case of (ii), a bath 
with a finitely many levels is sufficient while for (i), B must necessarily have infinite levels. 
The notation below follow the convention used in the proof of Theorem [31 

Sufficiency of (ii): Let Q be an operator that satisfies the conditions in (ii). Since ct/jq is 
isospectral with pQ, there is a unitary operator W on A such that 



For constructing the desired realization, the bath B will be chosen to be a system identical 
with A and its equilibrium state will be chosen to be pq. The realization is obtained by first 
transforming the device by W and then swapping the states of A and B. In other words, 
the reahzation isometry is Uab = S(W 1b) where the swap operator is defined to satisfy 
S \a)^ ® |«')^ = ® |«)^ for all |a) and \a'). If the initial state of the device is p, then 
the initial and final states of AB are 



(3Q = -W^{\npo + Sipo)tA)W 



(Bl) 



Pab = {p)a^Po , 

p'ab = {Po)a^{WpW^)b. 



(B2) 



(B3) 
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As a result, £{p) = po and the energy of the bath has been changed by 



AEb = ti Hb{WpW^ - Po) (B4) 

= -kBTtilnpoiWpW^ - Po) (B5) 

= -kBTtTp{WniipoW + S{po)lB) (B6) 

= tr pg (B7) 



and therefore the HTO associated with this reahzation is Q. This completes the proof of 
sufficiency of the special case (ii). 

Sufficiency of (i): Let Q be a hermitian operator satisfying only the strict inequality in 
the statement (i). The construction of the realization is quite similar to the construction 
used in the proof of Theorem [3l but now the bath has more components. 

Let d be the dimension of the Hilbert space of the device A and let r be the matrix rank 

n 

of po.[25]. Let X be a c?- level system and Y be an r- level system. The bath is taken to be 
composed of infinitely many copies of X and Y. The constitution of B can be expressed as 
B = (Y1Y2Y3 ■ ■ ■ )(XiX2X3 ■ ■ ■ ). Let the standard orthonormal bases of these two types of 
subsystems be denoted by the common notation where i = 0, 1, . . . , c? — 1 for X and 
i = 0, l,...,r — 1 for Y subsystems. The following shorthand will be used for the state 
vectors of the bath 

\J1J2 ■ ■ ■ ; ^1^2 ■■■)b = ® \3^)y, ® • • ■ ) ® (Ki)xi ® 1^2) X, ® ■ ■ ■ ) • (B8) 

The Hilbert space T-Lb oi the bath is defined as the linear span of these states where the 
associated sequences {jk} and {ik} contain only finitely many nonzero elements. These 
sequences can be considered as the digits of base-r and base-d representation of two non- 
negative integers m and n respectively. Namely, m = {■ ■ ■ iz32ii)r and n = {■ ■ ■ i'ii2ii)d- 
Therefore, I-Lb has a countable basis and it is a separable Hilbert space. 

Let and be the orthonormal sets in the Hilbert space of the device A that 

appear in the spectral decompositions of Q and po as in Eq. (129|) and 

■r-l 

Po = (PO > Pi > ■ ■ • > Pr-l > 0) (B9) 

where the eigenvalues are ordered in the indicated way. The isometry Uab on the composite 
system is defined as follows 

Uab ® \j1j2j3 ■ ■ ■ ; ^i^2«3 ■■■)b = l^ii)^ ® IhJsji ■■■ ] ^kii ■■■)b > (BIO) 
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FIG. 2: The isometry used in the reahzation for the complete erasure in the proof of the Theorem 
[2j The initial state in the top row is changed to the final state in the bottom row by copying the 
content of each subsystem to the next one on the right. 

for any allowed values of input labels. Essentially, the isometry Uab shifts the contents of 
the device and the subsystems in a chain which is infinite in two directions. The movement 
of the information among the subsystems is depicted in Fig. |2l Specifically, 

(i) the content of is copied to Yfc„i for each k >2, 

(ii) the content of Yi is copied to A in such a way that the standard basis of Yi is carried 
onto the eigenvectors of po, 

(iii) the content of A is copied to Xi in such a way that the eigenvectors of Q is carried 
onto the standard basis of Xi, 

(iv) and the content of X^ is copied onto X^+i for all k. 

Note that Uab maps an orthonormal basis of AB to an orthonormal set and therefore it is 
an isometry. 

For the Hamiltonian, the subsystems are taken to be independent and the Hamiltonians 
of individual subsystems are described in terms of a continuous parameter. For this purpose, 
consider continuous functions ei(z/) (i = 0, 1, ■ ■ ■ ,d — 1) and ej^u) (j = 0, 1, ■ ■ ■ , r — 1) of a 
real parameter u G [0, 1], to be used for energies of X and Y type subsystems respectively. 
They are chosen to satisfy the following properties. 
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(i) The zeroth level has zero energy, eo(i^) = eo(^) — 0) ^■ 

(ii) The energies of all excited levels diverge to +00 as the corresponding parameter ap- 
proaches to 1, i.e., 

lim eAu) = lim e'(i/) = 00 (for any i,j > 1). (BH) 
I/— >i I/— >i ■' 

(iii) The energy spectrum ej(i/ = 0) is identical with the eigenvalues of Q up to a constant 
shift as ei(0) — Qi — qo for all i. 

(iv) The energy spectrum e^ (z/ = 0) is such that the equilibrium density matrix is isospec- 
tral with po, i-e., e^(0) = kBT\n{pQ/pj). 

Let {z/fc} and {i^'k} be two infinite, increasing sequences starting from and converging 
to 1. The energy levels of and Y^. will be taken as ej(i/fe) and e^(t'fe) respectively. In other 
words, in terms of 

d-l 

h{u) = J2^ii^)\i){i\ > (B12) 

i=0 

h'{u) = Y,4{u)\i){i\ , (B13) 

the Hamiltonians of individual systems can be expressed as Hx,^ — h^Uk) and Hy^. — h'{v'i^. 
The bath Hamiltonian Hb is the sum of the subsystem Hamiltonians and therefore the 
energy eigenvalue of the bath state \j1j2 • • • ; ^1^2 " " " )b is 

fc=i 

Let = exp(— /3/i(z/))/C(z/) and cr'^u) = exp(— /3/i'(z/))/(^'(z/) be parameter-dependent 
density matrices and ({u) and C'{^) be the corresponding partition functions. The shorthand 
notations — cr(i/;t) and cr^ = cr'(^fe) ^^^^ be used for the states of and respectively. 
The partition function of the bath is then given by 

^B=ncK)c'K) ■ (B15) 

k=l 

Depending only on the large k behavior of the sequences {vk} and {ly'f,}, this product may 
converge or diverge. 

30 



In the following expressions, the subscripts are used for indicating the subsystem a given 
density matrix specifies to. Suppose that the initial state of the device A is p. The initial 
and final density matrices of the combined system are given by the following self-explanatory 
expressions. 

PAB = ■■■ {02)y2 ® (^i)yi ® {p)a ® (o^i)xi ® (c^2)x2 ■ ■ • 
p'ab = ■■■ i^3)Y, 8) {02)y, ® (po)a ® {p)x, ® {(ri)x, ■ ■ ■ 

It is obvious that the reduced quantum operation on the device is S{p) = po, i.e., the desired 
complete erasure. The change in the average energy of the bath is given by 

AEB = ti{pQ)~kBTiJ + Sipo)) + A (B16) 

where J = J{/3Q) and 

oo 

A = kBT{J + S{po))-tTaiQ + Y,iT{ak-i-ak)Hx, 

k=2 

oo 

+ J]trK+i-a^)/fy, , (B17) 

k=l 

which is a constant independent of the device state p. 
When Eq. (HM is used, A can be expressed as 

oo 

A = fc^T^(5K||afc+i) + 5K+ill^D) , (B18) 

k=l 

which shows that A is a strictly positive number. It will be argued below that the sequences 
of subsystem parameters {uk} and {z/^} can be selected in such a way that A can be adjusted 
to be equal to any positive value and this can be achieved with a finite Zb- As in the proof of 
Theorem [31 this will be done by showing that A can be made arbitrarily large and arbitrarily 
small. 

(1) First, if the parameter z/2 of X2 tends to 1, the final energy of X2 approaches to infinity. 
In this limit A diverges to 00. (ii) In the limit where all increments Uk+i — ^k and z/^_,_^ — z/^ 
tend to 0, the series in Eq. (]B17p approach to integrals which can readily be evaluated as in 
Eq. dM} and 

00 

tr(a^+i - a',)Hy^ -kBTS{a'{Q)) . (B19) 

k=l 
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Therefore, in this hmit, A approaches to 0. This shows that A can be chosen as small as 
possible. It can also be shown that the same can be done with a finite Zb value. 

Since A depends continuously on vy. and z/^, it is then possible to adjust these parameter 
sequences such that A has the value 

A = A;sT(J+^(po)) , (B20) 

which is a strictly positive quantity by our initial assumption. In that case, I\Eb = tr pQ 
and hence Q is the HTO associated with this realization. This completes the proof of the 
Theorem |2l □ 
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